High temperature correlation functions: universality, extraction of exchange 
interactions, divergent correlation lengths and generalized Debye length scales 
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We derive a universal form for the correlation function of general n component systems in the 
limit of high temperatures or weak coupling. This enables the extraction of effective microscopic 
interactions from measured high temperature correlation functions. We find that in systems with 
long range interactions, there exist diverging correlation lengths with amplitudes that tend to zero in 
the high temperature limit. For general systems with disparate long range interactions, we introduce 
the notion of generalized Debye length (and time) scales and further relate it to the divergence of 
the largest correlation length in the high temperature (or weak coupling) limit. 
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I. INTRODUCTION 

The study of correlation functions in systems with 
multi-component fields is an extremely general prob- 
lem having incarnations that range from systems in con- 
densed matter physics to fundamental field theories. One 
of the cornerstones of field theories and the study of criti- 
cal phenomena is the recognition of the universality that 
underlies general systems. This enables a unified un- 
derstanding and potent tools of analysis 1 3] . Any sys- 
tem generally displays a disordered high temperature (T) 
fixed point. Most of the research to date focused on the 
behavior of myriad systems at and in the vicinity of var- 
ious finite temperature transitions. In this work, we will 
focus on high temperature behavior and illustrate that a 
simple form of two-point correlation functions is univer- 
sally exact for rather general systems. This will enable 
us to make several striking observations. In particular, 
we will demonstrate that in contrast to common intu- 
ition, general systems with long range interactions have 
a correlation length that increases monotonically with 
temperature as T — > oo. As they must, however, the cor- 
relations decay monotonically with temperature (as the 
corresponding amplitudes decay algebraically with tem- 
perature). There have been no earlier reports of diverg- 
ing correlation lengths at high temperature. A thermally 
increasing length-scale of a seemingly very different sort 
appears in plasmas The Debye length, the distance 
over which screening occurs in a plasma, diverges, at high 
temperature, as Ajj cx VT. We introduce the notion of a 
generalized Debye length associated with disparate long 
range interactions (including confining interactions) and 
show that such screening lengths are rather general. 

Many early works investigated the high temperature 
disordered phase via a high temperature scries expansion 
[j| with an eye towards systems with short range inter- 
actions. In this paper, we report on our universal result 
for the Fourier transformed correlation function for sys- 
tems with general pair interactions. As it must, for near- 
est neighbor interactions, our correlation function agrees 



with what is suggested by standard approximate methods 
(e.g., the Ornstein-Zernike (OZ) correlation function that 
may be derived by many approximate schemes 0] ) • Our 
work places such approximate results on a more rigorous 
footing and, perhaps most notably, enables us to go far 
beyond standard short range interactions to find rather 
surprising results. Our derivations will be done for spin 
and other general lattice systems with multi-component 
fields. However, as illustrated later, our results also per- 
tain to continuum theories. 



II. OUTLINE 



In Section ITTTl we introduce the systems we study (gen- 
eral multi-component spin systems on a lattice; later 
sections generalize our result to other arenas - fluids, 
Bose/Fermi systems and so on). Section HVl contains a 
derivation of our main result about the universal form of 
the correlation function in the high temperature limit. In 
Section |V1 we comment on how the correlation lengths 
in a system behave in the high temperature limit. Sec- 
tion lVII introduces a generalized Debye length. In Section 
VIII we present some generalizations of our result to sys- 
tems which are not covered in Section Mil Section IVIIII 
outlines standard approximate techniques used to obtain 
our result. We give our concluding remarks in Section 

IS 



In Appendix[Al we show how to obtain a full high tem- 
perature series expansion of the correlation function to 
arbitrary order. In Appendix [B] we relate the general- 
ized Debye length over which long range interactions are 
screened to the diverging correlation length present in 
the high temperature limit. 
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III. SYSTEMS OF STUDY 

We consider a translationally invariant system with the 
Hamiltonian 



(i) 



The sites x and y lie on a d-dimcnsional hyper-cubic lat- 
tice with N sites having unit lattice constant. The quan- 
tities {S(x)} portray n-component spins ("O(n) spins") 
or general fields where ^(a;)! 2 = n at all lattice sites 
x. The normalization is adopted from Q. The case of 
n = 1 corresponds to Ising spins, n = 2 to XY spins, 
and so on. We assume that whenever the interaction 
kernel V(x) has a long range component, that compo- 
nent (unless stated otherwise) will always have some fi- 
nite screening, however small. This restriction is imposed 
to avoid well-known difficulties in taking thermodynamic 
limits in long range systems. In what follows, v(k) and 
Si(k) are the Fourier transforms of V{\x — y\) and Si{x). 
With this, Eq.flU) reads H = j W T,k v (k)s(k) ■ s(-k), 
up to an innocuous constant. Throughout, we employ 
the Fourier transform convention of a(k) = A(x)e zk x 
(and A(x) = j r j:^a(k)e-^ s ). 



IV. THE UNIVERSAL FORM OF THE HIGH 
TEMPERATURE CORRELATION FUNCTIONS 

We now derive a universal form for the correlation 
function at high temperature. As in any other calcu- 
lation with Boltzmann weights, the high temperature 
limit is synonymous with weak coupling. Initially, we 
follow standard procedures and examine a continuous 
but exact dual theory. High T (or weak coupling) in 
the original theory corresponds to strong coupling in the 
dual theory. We will then proceed to examine the conse- 
quences of the dual theory at high temperature where the 
strong coupling interaction term dominates over other 
non- universal terms that depend, e.g., on the number of 
components in the original theory. This enables an anal- 
ysis with general results. Unlike most treatments that 
focus on the character of various phases and intervening 
transitions, our interest here is strictly in the high temper- 
ature limit of the correlation functions in rather general 
theories of Eq.([T]). Our aims arc (i) to make conclusions 
concerning systems with long range interactions rigorous 
and (ii) to extract microscopic interactions from mea- 
surements. It is notable that due to convergence time 
constraints many numerical approaches, e.g., [||, com- 
pare candidate potentials with experimental data at high 
temperature (above the melting temperatures) where the 
approach that we will outline is best suited. We will per- 
form a transformation to a continuous but exact dual 
theory where the high temperature character of the orig- 
inal theory can be directly examined. 



We augment the right hand side of Eq.lQJ 
[— ^2g h(x) ■ S(x)] and differentiate in the limit h - 
to obtain correlation functions in the usual way. 

G(S-y) = -(S(x)-S(y)] 



by 
> 



lim 

h->o n 



1 



S 2 Z 



(2) 



with Z the partition function in the presence of the exter- 
nal field h. By spin normalization, G(x) = 1 for x = 0. 
The index i = 1,2, ...,n labels the n internal spin 
(or field) components. The partition function Z = 

Tr s [cxp (-^Ek<kMk)\2 + (3J2 s h(x) ■ S(x)J . 

The subscript S denotes the trace with respect to the 
spins. Using the Hubbard-Stratonovich (HS) transfor- 
mation, H [HI we introduce the dual variables {f]{x)} 
and rewrite the partition function as 



Z = Tr, 



/2 



drj^tye 2 ^ 



k(*0 \ +m(k) Si (-k) 



n 



fih{x)-S(x) 



(3) 



= A/"Tr. 



+ n ]r m ■ s(x) m ■ s(x) 



with U _1 (x) the inverse Fourier transform of l/v(k) and 
M a numerical prcfactor. The physical motivation in per- 
forming the duality to the HS variables is that we wish to 
retain the exact character of the theory (i.e., the exact 
form of the interactions and the 0(n) constraints con- 
cerning the spin normalization at all lattice sites). It 
is for this reason that we do not resort to a continuum 
approximation (such as that of the canonical 4 theory 
that we will discuss for comparison later on) where nor- 
malization is not present. Another reason to choose to 
work in the dual space is the correspondence with field 
theories that, in the dual space, becomes clearer in the 
high temperature limit (in which the quartic term of the 
^theories becomes irrelevant). Further details are in 
[ill- For O(n) spins, 



Z = N' 



Mr, 



exp 



n 



N 2 

^Y^V'^x-^rfix) ■ f)(y) 



n 



I n/ 2-i(V^\Nr}(x) + (3h(x)\) 
(y/n\NfK&) + /3K(x)|)"/2-i 



(5) 



The second factor in Eq. ([S]) originates from the trace over 
the spins and as such embodies the 0(n) constraints (the 
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trace in Eq.© is performed over all configurations with 
(S(x)) 2 = n at all sites x). Here, I u (x) is the modified 
Bessel function of the first kind. In the Ising (n = 1) case, 
the argument of the product in Eq.© is a hyperbolic 
cosine. Up to an innocuous additive constant, Eq.© 
corresponds to the dual Hamiltonian, 



H, 



■—tv-h 



-Vlnf - 



- y)rf(x) ■ fj{y) 
/2 -i(V^\Nrf(x) + PHx)\y 



■/3/i(x)|)"/2"i 



.(6) 



Our interest is in the h — > limit. The first term in Eq. 
(|()]) is the same for all n. This term dominates, at low /3, 
over the (second) n dependent term. As we will see, this 
dominance will enable us to get universal results for all 
n. From Eq.©, and the identity 



d_ 

dx 



Iu+l(x) 



we find that 



G(3-$ = 5 s , g + (1 - / ^f)'„lf\ x 

\\v(x)\\v(y)\ 

I n/2 (N^i\rf(x)\)I n/2 (N^\fim) \ 



^/2- 



^N^MxWnn-iiNMmi)/ / (7) 



where the average ((.)d) performed with the weights 
cxp(— PHj). Now, here is a crucial idea regarding our 
exact dual forms. From Eq.©, at high temperature, the 
variables rji{k) strictly have sharply peaked Gaussian dis- 
tributions of variance, 



-/3v(k) 



as /? 0. 



(8) 



Importantly, this variance tends to zero as /3 — > 0. By 
Parseval's theorem and translational invariance, 



-PV(Q)/N 2 



Thus, at high temperature, ((^(x)) 2 ) <C 1. It is therefore 
useful to perform a series expansion in the dual variables 
r\ and this would give rise to a high temperature series 
expansion in the correlation function. 

N 2 N 2 



N 1 



r]i{k) 



k.i 



X . i, 
2 



(9) 



with errors of C(l/T). Expanding Eq.© to C(l/T 2 ), 




(10) 



Eq.JTO]) leads to counter- intuitive consequences for sys- 
tems with long-range interactions. The second term in 
Eq. ([TU)) is independent of k and ensures that G(x) = 1 
for x = 0. Inverting this result enables us to find the mi- 
croscopic (spin exchange or other) interactions from the 
knowledge of the high temperature correlation function. 
We thus flesh out (and further generalize for multicompo- 
nent systems such as spins) the mathematical uniqueness 
theorem of Henderson for fluids [l2j for which a known 
correlation function G(x) leads to a known pair poten- 
tial function V{x) up to an innocuous constant. Eq. (|10[) 
leads to a correlation function which is independent of 
V(Q). Therefore, we can shift v(k) for all k's by an ar- 
bitrary constant or cquivalently set V(0) to an arbitrary 
constant. To 0(1/T), for V(0) = 0, we have, 



v(k) 



k R T 



k B T 



G(k) N G(fc') 



(11) 



The leading term of this expression for v(k) does not 
scale with T. This is because (1 — G(k)) oc 1/T at high 
temperature. Correlation functions obtained from exper- 
imental data can be plugged into the right hand side to 
obtain the effective pair potentials. Alternatively, in real 
space, for x ^= 0, 



V(x) 



-k B TG(x) + k B T G{x')G{x - x>) 



x'^Q,x 



(12) 



Note that the two terms in Ea.([T2"]) are 0(1) and 0(1/T) 
respectively, since G(x) is proportional to 1/T at high 
temperature for f ^ 0. Extension to higher orders may 
enable better comparison to experimental or numerical 
data. Our expansion is analytic in the high temperature 
phase (i.e., so long as no transitions are encountered as 
1 /T is increased from zero) . The Gaussian form of Eq. © 
similarly leads to the free energy density, 



2N ^ 



k B T 



v{k) 



+ 0(1/T). 



(13) 



Armed with Eq.(p?l [TU)) . we can compute any correlation 
function with the aid of Wick's theorem. For exam- 
ple, for unequal fc^s, we have, ((s(ki) ■ s(— fci))...(s(fc m ) • 

s(-k m ))} = (Nn) m Utx G (ki)- 

It is straightforward to carry out a full high temper- 
ature series expansion of the correlation function to ar- 
bitrary order. This is outlined in Appendix |XJ For ex- 
ample, to C(l/T 3 ), for V(x = 0) = 0, the correlation 
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function in real space, is given for i^O by. 
V(x) , 



G(x) 



k B T (k B T) 
1 



{k B Tf 



l —Y^V{z)V{x-z) 



2V(x) ^ V(z)V(-z) + 2 



n + 2 



.(14) 



V. HIGH TEMPERATURE CORRELATION 
LENGTHS 

We now illustrate that (i) in systems with short (or 
finite range) interactions, the correlation length tends to 
zero in the high temperature limit and (ii) in systems 
with long range interactions [l3j the high temperature 
correlation length tends to the screening length and di- 
verges in the absence of screening. 



A. Decaying lengthscales 

We consider first the standard case of short range 
interactions. On a hyper-cubic lattice in d spatial di- 
mensions, nearest neighbor interactions have the lattice 
Laplacian A(fc) = 2 J^ =1 (l — cosfc;) , with ki the l- 
th Cartesian component of the wave-vector k as their 
Fourier transform. In the continuum (small k) limit, 
A ~ \k\ 2 . Generally, in the continuum, arbitrary finite 
range interactions of spatial range p have v(k) ~ |fc| 2p 
with p > (and superposition of such terms thereof) 
as their Fourier transform. In general finite range inter- 
actions, similar multi-nomials in (1 — cosfc;) and in kf 
appear on the lattice and the continuum respectively. 
For simplicity, we consider v{k) ~ |fc| 2p . Correlation 
lengths are determined by the reciprocal of the imag- 
inary part of poles of Eq. (|T0)l . |Im We then 
have that in the complex k plane, (k*) 2p ~ —k B T. Poles 
are given by K ~ (/fc B T) 1 /( 2 P)exp[(2m+ l)m/(2p)\ with 
rn = 0, 1, 2p— 1. Correlation lengths then tend to zero 
in the high temperature limit as £ ~ J 1 ^ 1 /( 2 p) /| sin(2m + 
l)7r/(2p)| - there are p such correlation lengths. Simi- 
larly, there are p periodic modulation lengths scaling as 
L D ~ 27rT- 1 /( 2 rf/|cos(2m + l)n/{2p)\. The usual case 
of p = 1 corresponds to an infinite L B (i-e., spatially 
uniform (non-periodic) correlations) and £ ~ T~ 1//2 . 



B. Diverging lengthscales 

The novelty arises in the high temperature limit of sys- 
tems with long range interactions where v (k) diverges in 
the small k limit. Such a divergence enables the correla- 
tor of Eq. (p~0|) to have a pole at low k and consequently, 



on Fourier transforming to real space, to have a divergent 
correlation length. In the presence, of screening, v(k) di- 
verges and G(k) has a pole when the imaginary part of 
k is equal to the reciprocal of the screening length. The 
correlation length then tends to the screening length at 
high temperature. For concrctcness, we consider generic 
screened interactions where the Fourier transformed in- 

— with p 1 > and A 



teraction kernel vAk) ~ 775— 

the screening length. Perusing the poles of Eq.(fT0|). we 
find that for all p\ the correlation lengths tend to the 
screening length in the high temperature limit, 



lim e(T) 



X. 



(15) 



From Eq. (|15[) , when A becomes arbitrarily large, the cor- 
relation length diverges. Physically, such correlations en- 
able global "charge neutrality" [14[ for the corresponding 
long range interactions (Coulomb or other). This gen- 
eral divergence of high temperature correlation lengths 
in systems with long range interactions is related to the 
effective range of the interactions. At high temperature, 
the correlation function matches the "direct" contribu- 
tion, e -P y 'ff^ - 1 f3V eff (r). If the effective inter- 
actions between two fields have a range A, then that is 
reflected in the correlation length. In Coulomb systems, 
the Debye length, A^ sets the range of the interactions 
(for large distances, the interactions are screened). As 
stated earlier, at high temperature, A_d diverges. As seen 
by Fourier transforming Eq. (|10[) . though the imaginary 
part of the poles tends to zero (and thus the correlation 
lengths diverge), the prcfactor multiplying e™' 1 '^ is a 
monotonically decaying function of T. Thus in the high 
temperature limit the real space correlator G{x) mono- 
tonically decays with temperature (as it must). For in- 
stance, for p' = 1 in d = 3 dimensions, the pair cor- 
relator G(x) - e' x / x /{Tx) tends, for any non-zero x, 
to zero as T — > 00. [la ] That is, the amplitude van- 
ishes in the high temperature limit as (1/T). We find 
similar results when we have more than one interac- 
tion. For instance, in the presence of both a short and a 
long range interaction, (at least) two correlation lengths 
are found. One correlation length (or, generally, set of 
correlation lengths) tends to zero in the high tempera- 
ture limit (as for systems with short range interactions) 
while the other correlation length (or such set) tends to 
the screening length (as we find for systems with long 
range interactions). An example of a system where this 
can be observed is the screened "Coulomb Frustrated 
Ferromagnet" , fl5l - [l7| given by the Hamiltonian H = 
\-JY, {m S{x)S{y)+QY.^yVL(\x-y\)S(x)S(y% with 

J, Q > and the long range interaction Vl (x) = - — — 
in d = 3 dimensions and Vl[x) = Kq(x/X) in d = 2 with 
A the screening length and Kq a modified Bessel func- 
tion of the second kind. Similar dipolar systems [l8l - [2Cl | 
were considered. Apart from the usual correlation length 
that vanishes in the high temperature limit, we find an 
additional correlation length that tends to the screening 
length A. 
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VI. GENERALIZED DEBYE LENGTH (AND 
TIME) SCALES 

We now introduce the notion of generalized Debye 
length (and time) scales that are applicable to general 
systems with effective or exact long range interactions. 
These extend the notion of a Debye length from Coulomb 
type system where it is was first found. If the Fourier 
space interaction kernel v(k) in a system with long range 
interactions is such that is analytic at k = 0, then the 
system has a diverging correlation length, £i on g at high 
temperature. To get the characteristic diverging length- 
scales, we consider the self-consistent small k solutions 
to ksT /v(k) = —1 for high temperature (which gives the 
poles in the correlation function). Thus, as T — > oo, £i ong 
diverges as ^ksT, where p is the order of the first non- 
zero term in the Taylor series expansion of around 
k = 0. This divergent length-scale could be called the 
generalized Debye length. If the long-range interactions 
in the system are of Coulomb type, then this corresponds 
to the usual Debye length Xfj where p = 2. A more com- 
mon way to obtain this result is as follows. Suppose we 
have our translationally invariant system which interacts 
via pairwise couplings as in Eq.([l]). We can define a po- 
tential function for this system as, 



J2 V(\x-y\)S(y). 



(16) 



The "charge" S(x) in the system is perturbed by an 
amount S(x) and we observe the response <p{x) in the 
potential function <f>(x) assuming that we stay within the 
regime of linear response. We assume S(x) follows a 
Boltzmann distribution, i.e., S(x) = Aexp (— /3C<fi(x)), 
where C is a constant depending on the system. It fol- 
lows that S(x) = —f3CS(x)(j)(x). At this point, we can 
ignore the fluctuations in S(x) as it does not contribute 
to the leading order term. Thus, S(x) = —{3CSq4>{x), 
where So = (S(x)). In Fourier space, this leads to the 
relation, 

4>{k) = -pCS v(k)4>(k) (17) 
The modes with non-zero response are therefore given by, 



-v(k) oc ksT. 



(18) 



For a Coulomb system, these modes are given by 
(— fc~ 2 ) oc fcsT, yielding a correlation length Ajj oc 
y/k^T. 

As a brief aside, we remark that, repeating all of the 
above considerations (and also those to be detailed anew 
in Section lVlII A[) . if an imaginary time action for a com- 
plex field ip has the form 



- / dTdT'd d xd d x' 



ip(x,r) 

K{x - x' ,t - r')i>{x' ,t') \ + ..., 



(19) 



with the imaginary time coordinates < r, r' < (3 with 
a kernel K that is long range in space or imaginary time 
and the ellipsis denoting higher order terms (e.g., generic 
\tp\ A type terms) or imposing additional constraints on 
the fields ip (such as normalization that we have applied 
thus far for 0(n) systems) then the associated Debye 
length (or imaginary time) scale may diverge in the weak 
coupling (i.e., K — > aK with a — > + ) limit. In analogous 
way, repeating all of the earlier calculations done thus far 
for spatial correlations, we find that divergent correlation 
times in the low coupling limit for systems with a kernel 
K that is long range in |t — t'\. An action such as that of 
Eq. (fl9]) may also describe a system at the zero temper- 
ature limit (whence (3 — > oo) and the (imaginary) time 
scale is unbounded. 

In Appendix [B] wc will relate the divergence of the 
generalized Debye type length scales in the high tempera- 
ture limit to a similar divergence in the largest correlation 
length in systems with long range interactions. 



Confining potentials 

We discussed long range interactions (with, in general, 
a screening which may be set to be arbitrarily small) such 
as those that arise in plasma, dipolar systems, and other 
systems in condensed matter physics. In all of these sys- 
tems, the long range potentials dropped monotonically 
with increasing distance. Formally, we may consider gen- 
eralizations which further encompass confining potentials 
such as those that capture the effective confining po- 
tentials in between quarks in quantum chromodynamics 
(QCD) as well as those between charges in one dimen- 
sional Coulomb systems (where the effective potentials 
associated with the electric flux tubes in one dimension 
lead to linear potentials). The derivations that we car- 
ried throughout also hold in such cases. For instance, in 
a one-dimensional Coulomb system, the associated lin- 
ear potential V(x) ~ |x| leads to the usual Coulomb 
Fourier space kernel v(k) ~ k~ 2 . In general, for a po- 
tential V(x) ~ \x\~ a in d spatial dimensions, the corre- 
sponding Fourier space kernel is, as in the earlier case, 
v(k) ~ where p = d — a. Following the earlier dis- 

cussion, this leads, at asymptotically high temperatures 
(and for infinitesimal screening), to correlation lengths 
that scale as £ ~ \T. In the presence of screening, the 
correlation length at infinite temperature saturates and 
is equal to the screening length. Similarly, as seen by 
Eq. ([TB)) . the generalized Debye sreening length scales in 
precisely the same manner. In Eq. (|B2l) . wc will comment 
on the relation between the two scales. 



VII. GENERALIZATIONS 

Here we illustrate how our results can be generalized 
to systems which do not fall in the class of systems in- 
troduced in Section Hill 
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Disorder 



D. Bose/Fermi gases 



When Eq.([T]) is replaced by a system with non- 
translationally invariant exchange couplings V(x, y) = 
then V will be diagonal in an ortho normal ba- 
sis (| it)) different from the momentum space eigenstates, 
i.e., V\u) = v(u)\u). Our derivation will be identical in 
the \u) basis. In particular, Eq.([TU]) will be the same with 
v(k) replaced by v(u). 



B. Fluids 

Our results can be directly applied to fluids. In this 
case the spin at each site in Eq.((T]) may be replaced by 
the local mass density. The pair structure factor S(k) is 
the same as the Fourier space correlation function G(k) 



21 1. For r ^ 0, the pair distribution function g(r) is 



related to the correlation function G(r) defined above as 



g(r) = G{r) + 1. 



(20) 



For r = 0, g{r) = 0. 



C. General Multi-component Interactions 

In case of systems with multiple interacting degrees 
of freedom at each lattice site, we have a similar result. 
We consider, for instance, the non-rotationally invariant 
0(n) Hamiltonian, 



H = 7;^2^2 v ab(x,y)S a (x)S b (y), 



(21) 



x^y a, b 



where the interactions V a b (x, y) depend on the spin com- 
ponents 1 < a, b < n as well as the locations x and y. 
By fiat, in Eq.([2"Tj). V a b(x = y) = 0. Non-rotationally 
symmetric interactions such as those of Eq. (|2"Tj) with a 
kernel V a b which is not proportional to the identity ma- 
trix in the internal spin space 1 < a, b < n appear in, e.g., 
Dzyaloshinsky-Moriya interactions [22| , isotropic 12311 and 
non- isotropic compass (24[ , Kugel-Khomskii [2^, [23] and 
Kitaev type [26j models. Such interactions also appear 
in continuous and discretized non-abelian gauge back- 
grounds (and scalar products associated with metrics of 
curved surfaces) used to describe metallic glasses and 
cholesteric systems [27rf3l| . The lattice "soccer ball" spin 
model [27| is precisely of the form of Eq. (f2"Tj) . Replicat- 
ing the calculations leading to Eq. (|A6[) . for x ^ y, to 
C(l/T 2 ), we find that 



Gab{x,y) = (S a (x)S b (y)} 
V a b(x,y) 1 
k B T {k B Tf 



Y,V ac (x,z)V cb (z,y).(22) 



Here we discuss Bose/Fermi systems to illustrate the 
generality of our result from Eq. (fT0| We consider the 
Hamiltonian given by 



H = H + Hi, 



where 



X 



(23) 



with p(x) = ip* (x)'ip(x) — (V'^V')o- 

Here and throughout, (-)o denotes an average with re- 
spect to Ho (the ideal gas Hamiltonian). The fields ip 
obey appropriate statistics (Bose-Einstein/Fermi-Dirac) 
depending on the system being studied. The standard 
partition function is 



Z = Zq Dr)(x, t) e 



(24) 



Here, r is the standard imaginary time coordinate (0 < 
T < 0) ■ Zqis the partition function of the non- interacting 
system described by Hq, 77-s are the dual fields after per- 
forming the HS transformation. We can express $ as 



iV 2 f 13 

x,x' 

~In^T r exp^ijT / 'dT^» 7 (f,T)p(£,T)^ . (25) 

where T T is the (imaginary) time-ordering operator. It 
is clear that the factor of the partition function which 
controls high temperature behavior comes from the first 
term in Thus, for small (3 (high temperature), the 
distribution of the values of 77 is sharply peaked around 
zero. Also, for small /3, the integrands of Eq. (|2l)]) have 
little dependence on r. Therefore, at high temperature, 



N 7 



2^ 2 ^ 



£>(a0 [V~\x-x') 



+PA(x-x , )]t } {&), 



(26) 



where A(x — x') = (p(x)p(x')) a = CSggi, with C = p\ 
being a constant. The correlation function for this system 
is defined as Gix — y) = (p(x)p(y)), It is easy to show 
that written in terms of the dual variables, 



G(x-y) 



/ f'(N V (x))f(N V (y)) \ 
\f(NTj(x)) f(N V (y))/ d 



(27) 



where /(a) = Tr^^ e ap ^ and, as before, (-)d denotes 
the average with respect to the dual fields rj. For small 
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values of the rj variables (high temperature), we have in 
general, G(x — y) = C + Ci(rj(x)ri(y))d, with C chosen 
such that G(x) — C for x = and C\ defined by the 
statistics of p and the form of the pair interaction V. 
Therefore, we have, 



G{k) = C + 



C[Cv(k) + k B T] 

v 2.^ 



N Y C[Cv(k) + k B T] 



(28) 



This is similar to the classical 0(n) correlation func- 
tion in Fourier space [Eq. (fTU|) ] We can easily generalize 
Eq. (|28|) for multi-componcnt/polyatomic systems as in 
Eq. (|22[) . Applied to scattering data from such systems, 
our results may enable the determination of effective un- 
known microscopic interactions that underlie the system. 
Similarly, replicating the same derivation, mutatis mu- 
tandis, for quantum SU(2) spins S = (S x , S y , S z ) in the 
coherent spin representation leads to the high tempera- 
ture result of three-component (0(n = 3)) classical spins. 
This illustrates the well known maxim that at high tem- 
perature, details may become irrelevant and systems "be- 
come classical" . In a similar manner, at high T, the de- 
tails underlying the classical 0(n) model (the 0{n) nor- 
malization constraints concerning a fixed value of 
for n component vectors S(x) at all sites x) effectively 
became irrelevant at high temperature - the behavior for 
all n was similar. 



A. Ginzburg-Landau 4 -type theories 

In the canonical case, the free energy density of the <fi 4 
theory is given by 

= + \r^(x) + ^(x). (29) 

A finite value of a corresponds to the "soft-spin" approx- 
imation where the norm is in not constrained, (4> 2 (x)} ^ 
1. Here, r = c(T — Tp), with c a positive constant. 
The partition function [38[ is Z = j D(f> e~ F where 
F = j J 7 d d x with d the spatial dimension. At high 
temperature, the correlator behaves in a standard way 
(the OZ form) (J0(fc)| 2 ^ = -p^. The irrelevance of 

the </> 4 term may, e.g., be seen by effectively setting 
4> 4 {x) — > Q((j) 2 (x))(f) 2 (x) in the computation of the par- 
tition function. As (4> 2 {x)) is small [in fact, from Fourier 
transforming the above, ((j) 2 (x)} = 0(1/T)], the c/> 4 term 
is smaller than the (V<t>) 2 term in Eq. (f2"9"|) by a factor 
of a/T and therefore can be neglected. When general 
two body interactions with an interaction kernel v(k) are 
present, wc similarly have (\<fi(k)\ 2 ) = J + . Our result 

of Eq. ([10p for interactions of arbitrary spatial range is 
new and illustrates that suggestive results for the corre- 
lation lengths attained by soft spin approximations are 
not far off the mark for general systems in the high tem- 
perature limit. As far as we are aware, the high tempera- 
ture correlation length of general theories was not known 
to be similar to that suggested by various perturbative 
schemes (including the 1/n [3!| and e expansions (40|). 



B. Correlation Functions in the large n limit 



VIII. APPROXIMATE METHODS 



The exact high temperature results that we obtained 
for lattice spin systems and the generalizations that we 
discussed in Section [VIII arc, as wc will show below, sim- 
ilar to those attained by several approximate methods. 
This coincidence of our exact results with the more stan- 
dard and intuitive approximations enables a better un- 
derstanding from different approaches. A corollary of 
what we discuss below is that the divergence of the cor- 
relation lengths in systems with long range interactions 
in the high temperature limit (as in Section |Vl) appears 
in all of these standard approximations. However, as we 
illustrated earlier in our work, and in Section [V] in par- 
ticular, this divergence is not a consequence of a certain 
approximation but is an exact feature of all of these sys- 
tems in their high temperature limit. 

In what follows, we will specifically discuss (i) </> 4 field 
theories, (ii) the large n limit, and (iii) the OZ ap- 
proach for fluids invoking the mean-spherical approxi- 
mation (msa) [13. 



We now provide a derivation of Eq. (jTD| as it applies 
in the large n limit. Long ago, Stanley [7| demonstrated 
that the large n limit of the 0(n) spins is identical to the 
spherical model first introduced by Berlin and Kac |4lj . 

The single component spherical model is given by the 
Hamiltonian, 



H 



x-y\)S(x)S(y). 



(30) 



x=£y 



The spins in Eq.([30| satisfy a single global ("spherical") 
constraint, 



Y,S 2 (x) = N, 



(31) 



enforced in its average value [37J by a Lagrange multiplier 
/x. This leads to the functional H' = H + fj,N which 
renders the model quadratic (as both Eqs. (f5Ul f3TT) are 
quadratic) and thus exactly solvable, see, e.g., f fl6l|). 

From the equipartition theorem, for T > T c , where no 
condensate is present, the Fourier space correlator 



1 



G(k) = -(\s(k)\ 



k B T 
v{k) + fx 



(32) 
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The real space two point correlator is given by 



IX. CONCLUSIONS 



N 



(33) 



To complete the characterization of the correlation func- 
tions at different temperatures, we note that the La- 
grange multiplier /z(T) is given by the implicit equation 
1 = G(x = 0). Thus, 



E 



i 



N ^ v (k) + fj, 



(34) 



This implies that the temperature T is a monotonic in- 
creasing function of p. Eq. (|3"4")) also implies that in the 
high temperature limit, 



knT. 



(35) 



Taken together, Eqs. (|32l35p yield Eq.flTU]) in the asymp- 
totic high temperature limit. For completeness, we 
briefly note what happens at low T (T < Tc)- In the 
spherical model, at the critical temperature (T c ), the La- 
grange multiplier p takes the value, 



min{u(fc)}. 

k 



(36) 



For T < T c , (at least) one mode q is macroscopi- 
cally occupied, the mode(s), q being occupied is one 
for which v(k) is minimum. The "condensate fraction" 
(| S (<z)| 2 )/iV 2 >0. 



C. Ornstein-Zernike Equation 

As noted earlier, application of the MSA to the OZ 
equation for fluids reproduces similar results for the 
"total correlation function", h(f). This is defined as 
h(r) — g{f) — 1, where g(f) is the standard radial dis- 
tribution function. The OZ equation for a fluid with 
particle density p is given by, 



h(f) = C(f) + p / dr'C{r- ?)h{?), 



(37) 



where C{r) is the "direct correlation function". Using 
the MSA, C(r) = — 0V(r) [37|, we get in Fourier space, 



S(k) 



k R T 



pv(k) + ksT 



(38) 



This is similar to our result for G(k). However, it is valid 
only for systems in which the MSA is a good approxima- 
tion. 



(i) We derived a universal form for high temperature 
correlators in general 0{n) theories. This enables the ex- 
traction of unknown microscopic interactions from mea- 
surements of high temperature correlation function. 

(ii) We discovered divergent correlation lengths in sys- 
tems with long range interactions in the high temperature 
limit. This divergence is replaced by a saturation when 
the long range interactions are screened. 

(iii) Wc introduced generalized Debye lengths associ- 
ated with such divergent correlation lengths. 



Appendix A: High temperature series expansion of 
the correlation function 



We now outline in detail how we may obtain a high 
temperature (T) series expansion of the correlation func- 
tion to arbitrary order for a general system with transla- 
tional invariance. The result provided above was derived 
to order C(l/T 2 ). This and the results we present below 
are valid in the high temperature phase of general lattice 
(spin or other) and continuum systems. However, it will 
hold in lower temperature phases of the system provided 
we can analytically continue to those phases from the 
high temperature phase, i.e. attain those phases without 
having a phase transition. It is also worth mentioning 
that since we set the temperature to be arbitrarily high, 
the density does not have to be small as is assumed in 
methods derived from Mayer's cluster expansion for flu- 
ids. Our result is therefore valid for the high temperature 
phase of any system. In general, the long range charac- 
ter of the interactions will not enable us to invoke many 
of the simplifying elegant tricks present elsewhere. For 
instance, the counting of connected contours and loops 
[42l [43| that appear in high temperature series expansion 
involving nearest neighbor interactions cannot be applied 
here. 

We can perform the high temperature series expansion 
directly in the original spin space. However, we find it 
easier to make a transformation to a dual space where our 
Boltzmann weights become Gaussian in the high temper- 
ature limit. 

The correlation function of the original theory can be 
expressed in terms of the correlation function (and higher 
moments) of the dual theory - we employ that in our 
calculation. The dual theory to a nearest neighbor fer- 
romagnetic system is a Coulomb gas. Nearest neighbor 
ferromagnetic system in dimensions d > 2 at low T has 
an ordered phase and a small correlation length (corre- 
lation length diverges at T = T c ). This does not imply 
that the Coulomb system has a small correlation length 
at high temperature. 0(n) constraints become faint at 
high T in the dual theory whereas in the exact Coulomb 
gas at high T, the 0(n) constraints are there. The same 
also applies for a soft spin realization of the Coulomb gas 
where exp[— /3u(S 2 — l) 2 ] which is zero as (3 — > oo (or 



9 



T —> 0) unless S 2 = 1 everywhere. By contrast in the ex- 
act dual theory at high temperature, the relative strength 
of the 0(n) constraints becomes negligible relative com- 
pared to the "interaction" term containing (0V)~ 1 . Even 
though we can ignore prefactors when = 0(1) and 
consider dual theories and soft spin realization we cannot 
ignore the T dependence at high T about the infinite T 
disordered limit. Otherwise we get a contradiction as our 
exact calculation with the exact dual theory (containing 
the T dependent prefactors) shows. 

We will keep things general and perform the simple 
series expansion of the dual Hamiltonian Hd in Eq. (J5J) ■ 

N 2 



0^ 



4/2-i(V^v|#e)I) 

(^iN\ff(x)\) n / 2 - 1 J ' 



20 2 
N 2 



N 



A(n + 2)0 



Y ■ rf(x)Y 



(Al) 



In Eq. (|Aip . the interaction V should be thought of as 
a translationally invariant matrix. That is, in a Dirac 
type notation, (x\V\y) = V(x — y). In Eq. (|Al|) . V~ x is 
the inverse Fourier transform of l/v(k), where v(k) is the 
Fourier transform of V(x) (44[. 

Next, we separate Hd into a quadratic part Hdo and 
higher order (interaction type) terms which we denote by 
AH. That is, 



H, 



(10 



N 2 



AH 



X 1 



4(n + 2)0 ^ 



(A2) 
(A3) 



The expectation value of any quantity X may be com- 
puted by 



(X) d = 



-/3AH> 



(A4) 



= (X) do -f3[(XAH)do-(X) do (AH) do ] 
B 2 

+ |y [(X(AH) 2 ) d() - 2(XAH)do(AH) do 

+2(X)do{AH) 2 - (X) d0 {(AH) 2 ) d0 ] 

+... , (A5) 

where (-)do represents the expectation value calculated 
with the Boltzmann weight associated with the Hamilto- 
nian H d Q. We may retain terms to arbitrary order in r/ 2 



(or corresponding order in 1/T). Eq.([7]) can be expanded 
to arbitrary order in r] 2 where we rewrite all expectation 
values with respect to the Hamiltonian Hdo- The terms 
become expectation values of a product of an even num- 
ber of r\ fields with respect to the quadratic Hamiltonian 
Hdo- We can then use the Wick's theorem to compute 
the expectations with respect to Hdo to all orders. To 
order 1/T 3 we obtain for x ^ 0, 



G{x) 
1 



{k B Tf 



k B T 



2V(0)V(x) 



J2v(y)V(z)V(x-y-z) + 
2V(x) V(z)V(-z) + 3V(0) ^ V(z)V(x - z) 

z z 

(V{x)f 



-5(V{0)yV(x)-2- 



n + 2 



(A6) 



As a brief aside, we note that from the fluctuation dis- 
sipation theorem the susceptibility x = 0TlisG{x). At 
asymptotically high temperature, G(x) ~ 5g t o giving rise 
to Curie's law, \ cx 1/T. The terms in Eq. (|A6j) lead to 
higher order corrections. To next order, x = k B (T-e c ) 
with the Curie temperature 9c = J2s^o V(%) m * ne wea k 
coupling limit. Thus far, in the literature, the Curie- 
Weiss form was invoked to ascertain whether a given sys- 
tem has dominantly ferromagnetic or anti-ferromagnetic 
interactions (sign of 9c) and their strength (|#c|)- We 
see that by not focusing solely on \ ~ 0G(k — 0) but 
rather on the scattering function G(k) for all fc, we can 
in principle deduce the interaction v(k) and hence V(x). 

In Fourier space, the real space convolutions become 
momentum space products and vice versa. Eq. (|A6[) then 
reads 



G{k) = 1 - 



+ 



{k B Tf 



(k B T) 

Mk)? + 2 -^ 



fel 



-3V(x = 0)(v(k)) 2 - 5(V(x = 0)) 2 v(k) 



N 2 (n + 2) 



^ v(ki)v(k 2 )v(k - ki - k 2 ) 



k\ ,k-2 



Gi(0), 



(A7) 



where Gi(0) is the value obtained by inserting x = 
in Eq. (|A6[) . It should be noted that the real space cor- 
relation function cannot change if we shift the on-site 
interaction V(x = 0) which is equivalent to a uniform 
shift to v(k) for all k. This is because the 0(n) spin is 
normalized - jS^ic)! 2 ~ n at all sites x. This invariance 
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to a constant shift holds for all T and consequently to 
any order in 1/T, the coefficients must be invariant to 
a global shift in v(k). Amongst other things, we earlier 
invoked this invariance [ll| to shift v(k) to enable a HS 
transformation in the cases for which initially v(k) > 
for some values of k. We can, of course, invoke this in- 
variance also here to obtain the above high temperature 
series expansion with a well defined HS dual. The final 
results, as we re-iterated above are invariant under this 
shift as is also manifest in our series expansion in powers 
of 1/T. Although obvious, we note that the expansion 
in Eq. (|A7[) is performed in power of 1/T involving v(k) 
for real vectors k. In examining the correlation lengths 
via contour integration in the complex k plane, the cor- 
responding v(k) may be extended for complex k. 

We see from the expansion in Eq. (|A7[) that already to 
0(1/T), it is also clear that the lengthscales of the system 
(which are determined by the poles of the Fourier space 
correlation function) are governed by the poles of v(k) in 
the complex k space. Thus, if, e.g., v(k) = l/(k 2 + A~ 2 ), 
the correlation length tends to A at high temperature. It 
therefore must diverge for a system with no screening. 

In cases where the correlation function is known from 
some experimental technique or otherwise, the series ex- 
pansion for the correlation function can be inverted to 
arbitrary order to obtain the pairwise interactions. To 
0(1/T 2 ), for non-zero separation x, the potential func- 
tion is given by, 



V{x) = -k B T 



z 

J^G(y)G(z)G(x-y- 



y,z 



+2G{x)J2'g(z)G(-z) 



2(GQ?)r 
n + 2 



(A8) 



The prime indicates that the sum excludes terms con- 
taining G(0). As is evident from our earlier results and 
discussion, in Eq. (|A8p . each correlation function G(x) is 
of order (1/T). 

We re-iternate that as in our discussion in Section 
IVIII1 our results for lattice 0(ri) spin models match with 
the leading order behavior at high temperature obtained 
from several standard approximate theories based on 
Mayer's cluster expansion derived for liquid systems, 
e.g., Born-Green theory [45| and OZ theory with Percus- 



Yevick approximation [46( or MSA 37 1. As implicit 
above, our 1/T expansion can indeed be extended to 
systems in which the liquid and the gas phase are not 
separated by a phase transition, e.g., for pressures larger 
than the pressure at the liquid- vapor critical point. As 
further noted in Section IVIII1 various approximations 
also suggest that at high temperature, the correlation 
length may match the length-scale characterizing of the 
interaction potential and, in particular, would diverge 
in systems having long range interactions (as we have 
established). 



Appendix B: Relation between the generalized 
Debye lengths and divergence of the high 
temperature correlation lengths 

An intuitive approximate approach for the understand- 
ing of the rigorous yet seemingly paradoxical result that 
we report in this work- that of the divergence of the cor- 
relation lengths in the high temperature limit of systems 
with long range interactions is afforded by the OZ frame- 
work. Specifically, in the language of OZ approximations, 
the "total" high temperature correlation function is the 
same as the "direct" correlation function (see, e.g., plj 
(section 2.6) for the definition of the "direct" OZ corre- 
lation functions) and behaves as 



G{x) ~ -PV(x) 



(Bl) 



for x ^ 0. Thus, if the potential is screened beyond a 
distance A, the correlation length approaches A at high 
temperature. That is, if we have an effective interac- 
tion resulting, e.g., from higher order effects in 1/T, such 
as that leading to the Debye screening length (Xd) in 
Coulomb systems (and generalizations introduced earlier 
in Eg . (fT5|) ) . then at high temperature, the correlation 
length 



■+ A, 



(B2) 



This is a particular case of Eq. (|T5]) . 

To C(l/T 2 ), Eq.flM} is identical to Eq.^TDJ). The poles 
of G in the complex k plane can, of course, be computed 
to by finding those of Eq. (fTU)) or considering those di- 
rectly of Eq. (|A7[) : both give rise to the same answer as 
they must. 
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